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In recent years, unconventional metamaterial properties have triggered a revolution of electromag¬ 
netic research which has unveiled novel scenarios of wave-matter interaction. A very small dielectric 
permittivity is a leading example of such unusual features, since it produces an exotic static-like 
regime where the electromagnetic held is spatially slowly-varying over a physically large region. The 
so-called epsilon-near-zero metamaterials thus offer an ideal platform where to manipulate the inner 
details of the ’’stretched” held. Here we theoretically prove that a standard nonlinearity is able 
to operate such a manipulation to the point that even a thin slab produces a dramatic nonlinear 
pulse transformation, if the dielectric permittivity is very small within the held bandwidth. The 
predicted non-resonant releasing of full nonlinear coupling produced by the epsilon-near-zero con¬ 
dition does not resort to any held enhancement mechanisms and opens novel routes to exploiting 
matter nonlinearity for steering the radiation by means of ultra-compact structures. 


During the last decade the metamaterial route for 
achieving unusual electromagnetic properties has at¬ 
tracted a great deal of interest in both theoretical and 
applied research, since it has brought to light novel elec¬ 
tromagnetic regimes [H-Q, and has suggested a number 
of remarkable devices for extreme manipulation of the 
radiation [^-Q. Structures exhibiting very small dielec¬ 
tric permittivity, or epsilon-near-zero (ENZ) metamate¬ 
rials [m, belong to the family of media able to af¬ 
fect electromagnetic radiation in a very unconventional 
way because the medium’s effective wavelength is much 
larger than the vacuum wavelength, and because they 
host a regime where both field amplitude and phase 
are slowly-varying over relatively large portions of the 
bulk, which is quite opposite to geometrical optics. Such 
key feature has been exploited to conceive setups where 
ultra-narrow ENZ channels are able to ’’squeeze” elec¬ 
tromagnetic waves at will (l^ - [r^ . and to develop new 
paradigms of devices for tailoring the antenna radiation 
pattern [Tsl - I^ . In addition, ENZ metamaterials have 
also been shown to support a rich phenomenology of sur¬ 
face waves to achieve perfect absorption (2^, to 

enhance spatial dispersion effects [2^, and to support 
novel cloaking mechanisms SEl. 

Other interesting mechanisms and effects arise when 
the ENZ regime is combined with matter nonlinearity. A 
kind of crucial interplay between small linear permitti vity 
and optical nonlinearity has been identified in Ref. |32l| . 
where the authors discuss an all-optical transition from 
dielectric to metal behavior in a nonlinear ENZ slab. 
A similar mechanism has been exploited to predict a 
class of solitons [3^ where the intensity-driven dielectric- 
metal transition occurs along the transverse soliton pro¬ 


file, yielding exotic features like transverse power flow 
reversing [^, unusual shapes like hollow-core [s^, and 
two-peaked and flat-top profiles [s^. Eurthermore, the 
combination of the ENZ regime with nonlinearity ben¬ 
efits from the non-resonant enhancement of the normal 
electric field component across the vacuum-ENZ medium 
interface [s^, producing intriguing effects like transmis¬ 
sivity directional hysteresis j38l . l39| and enhancement of 
second and third harmonic generation (40l - l43| . A dif¬ 
ferent field enhancement mechanism has been identified 
within narrow ENZ plasmonic channels, which has been 
exploited to boost optical nonlinearities [i^, to inves¬ 
tigate temporal soliton excitation [d^, and for the en¬ 
hancement of second-harmonic generation efficiency (d^ . 

In this work we theoretically prove that the interplay 
between the ENZ condition and the optical nonlinearity 
triggers a novel nonlinear matter-wave coupling which 
literally unlocks the full potential of the generally weak 
matter nonlinear response. By theoretically investigat¬ 
ing the scattering of electromagnetic pulses by a thin 
nonlinear dispersive slab we show that a marked non¬ 
linear pulse dynamics occurs only if the absolute value 
of the dielectric permittivity is very small over the pulse 
bandwidth. Such a nonlinear scenario stems from the fact 
that, in the ENZ regime, the field is spatially slowly vary¬ 
ing and accordingly not characterized by a large number 
of ’’nodes” around which its amplitude is small, as the 
medium nonlinearity is allowed to uniformly affect the 
field over the entire bulk. In other words, as a viable 
strategy to attain a highly nonlinear response, here we 
suggest to enlarge the physical volume over which the 
nonlinearity is effective by means of the ENZ condition, 
as opposed to standard approaches that resort to field 


2 



FIG. 1: Pulse scattering by a slab supporting the ENZ regime, (a) Geometry of the scattering interaction between 
a transverse magnetic pulse which is both spatially and temporally localized and a material slab having, in the linear regime, 
zero-crossing-points of the real part of its dielectric permittivity, (b) In the linear regime where the medium polarization 
|P| is much smaller than the saturation polarization Ps (see Eq. O), the slab has a dielectric permittivity e(tu) (see Eq.®) 
with a standard Lorentz profile located at the resonant frequency oJe and with a zero-crossing-point of its real part at ujo (see 
Eq.Q). Dispersion parameters have been chosen in such a way that the imaginary part of the permittivity is low around 
zero-crossing-point so that |e(a;)| is much smaller than one in a spectral bandwidth around a;o and the slab can support the 
ENZ regime. 


enhancement mechanisms or giant nonlinear parameters. 


I. RESULTS 


Pulse scattering by a slab supporting the ENZ 
regime. The ENZ regime occurs when |e(a;)| <C 1 holds 
over the entire spectral bandwidth of the incident field. 
In conjunction with low absorption, this occurs if the field 
has a relatively narrow-band whose frequencies are very 
close to a zero of the real part of the medium permit¬ 
tivity. Permittivity zero-crossing-points are present both 
in metamaterials, where they can be tailored by choos¬ 
ing the relative metal-dielectric content of the unit cell, 
and in standard materials close to the resonant absorp¬ 
tion frequencies. Here we analyze the second situation, 
and we model the dynamics of the matter polarization 
P driven by the radiation electric field E through the 
equation 

P = eo (e. - 

( 1 ) 

where Ps is the saturation polarization that governs non¬ 
linear oscillator behavior. For |P| very much smaller 
than Ps, Eq. © reduces to the single-pole Lorentz os¬ 
cillator with resonant frequency cOg, loss coefficient SgUJe 
and static dielectric permittivity Cg. Therefore, in the 
linear regime, the dielectric permittivity experienced by 
monochromatic fields is 


e{u}) = 1-1- 



( 2 ) 


and it admits the zero-crossing-point 



(Cs -I- 1 — < 5 g ) — 4:€s 




( 3 ) 

since Re[e(a;o)] = 0. For larger |P| such that the con¬ 
dition |P| <C Ps still holds, Eq. © describes a stan¬ 
dard Kerr response whereas, at higher polarization 
strengths, it accounts for a realistic saturation of matter 
nonlinear response [s^ . Therefore, Eq. © contains all 
the ingredients necessary to describe a realistic and very 
general nonlinear wave-matter interaction in the ENZ 
regime. 

In order to investigate such a regime we consider the 
scattering interaction reported in Fig.la, where an elec¬ 
tromagnetic pulse is launched from vacuum to orthogo¬ 
nally impinge on the surface of a dielectric slab whose 
polarization is described by Eq. ©. The material dis¬ 
persion parameters are Sg = 0.01 and Cg = 1.2 so that 
uiQ = 1.095we and the imaginary part of the permittiv¬ 
ity around the zero-crossing-point is small, as reported 
in Fig.lb and in its inset, thus allowing the ENZ regime. 
We have set for the slab thickness L = 1.25Ae where 
Ae = 27rc/a;e is the resonant wavelength and there¬ 
fore, since we will mainly consider pulses whose band¬ 
width is localized around the zero-crossing-point, the 
slab thickness is comparable with the main field wave¬ 
length. The pulse is a transverse magnetic (TM) excita¬ 
tion whose profile at the launching plane Ex{x, Zin,t) = 


Eoe~'^e sin(u)t) has time shift to = 3.178T0^a;“^ 

and it is both spatially and temporally localized, a = 
1.25Ae and r = 1.059 • 10^a;“^ being its transverse and 
temporal widths, respectively. The amplitude Eq and 
the carrier frequency cD will be varied to investigate the 
ENZ regime. All the situations considered below are 
in the quasi-monochromatic regime since the carrier fre¬ 
quency uj is always comparable with Ug and the pulse 
spectral width i5a; ~ 1/r = 9.442- lO^'^We is much smaller 
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FIG. 2: Nonlinear signatures of the scattering in the ENZ regime, Gj = ujq. (a) Spectrum \e^\ of the incoming 
pulse (arbitrary units) and real and imaginary parts of the slab linear dielectric permittivity e. The interaction is in the ENZ 
regime since |e| <C 1 holds over the pulse spectrum, (b) Absolute values of the dimensionless analytic signals, normalized to 
the amplitude of the incoming pulse eg, |ef|/eo and |ef|/eo of the components of the outgoing pulse at the exit slab plane 
Zout = ciJeL/(2c). The dependence of |ef |/eo and |ef|/eo on eg is the key signature of the nonlinear wave-matter interaction, 
(c) Square absolute values of the dimensionless analytic signal of the polarization |pf -I- |pf P, dimensionless displacement field 
components spectra \d^\ and | and dimensionless electric field components spectra multiplied by permittivity absolute value 
\te^\ and |eef | at the slab middle plane Z = 0. The occurrence of the nonlinear wave-matter interaction is further testified 
by two facts: 1) -I- |pfP is not very much smaller than one and 2) \d^\ ^ |eef^| and |df | ^ |eef |, i.e. the standard linear 

constitutive relation does not hold. 


than We. Finite difference time domain simulations have 
been performed to solve Maxwell equations coupled with 
Eq. dl]) (see Methods). Results are reported using the di¬ 
mensionless coordinates iX,Z) = uje{x,z)lc, T = uj^t 
and fields ^z) / Ps-i 

{px,Pz) = {Px,Pz)/Ps which remove the resonant fre¬ 
quency We and the saturation polarization Pg from the 
model (see Methods). 

Highly nonlinear wave-matter interaction in the 
ENZ regime. In the first set of simulations we have con¬ 
sidered various pulses with carrier frequency located at 
the zero-crossing-point, i.e. w = wq, with different ampli¬ 
tudes eo = egEg/Pg. In Fig.2a we have plotted the abso¬ 
lute value of the Fourier transform (arbitrary units) 
of the incoming pulse ex{0, Zin,T) superimposed to the 
real and imaginary parts of the permittivity. This fig¬ 
ure clearly shows that |e(w)| is much smaller than one all 
over the field spectrum, thus proving that the interaction 


occurs in the ENZ regime. 

In Fig.2b, which contains the main result of the present 
paper, we report, for a number of different amplitudes eg, 
the space-time profiles of the outgoing pulses. Specif¬ 
ically we plot the absolute values of the analytic sig¬ 
nals ef(X, Zout,T) and ef{X,Zout,T) of the field com¬ 
ponents, normalized with respect to the amplitude eg, at 
the exit slab surface Zout = WeT/(2c). For eg = 10“® 
and eg = 10“^ the normalized outgoing pulses are iden¬ 
tical. This shows that the pulse amplitudes are small 
enough to prevent the matter nonlinearity from influenc¬ 
ing the scattering. By increasing the pulse amplitude, 
for eo = 2 • 10“^ and eo = 4 • 10“^, the normalized out¬ 
going pulses undergo dramatic transformation. The de¬ 
pendence on the input intensity is a key signature of the 
onset of the nonlinear regime. At higher pulse intensi¬ 
ties, the normalized outgoing pulses having amplitudes 
eo = 1.2 and eo = 1.6 are almost identical, so that the 
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FIG. 3: Linear character of the scattering outside of the ENZ regime, u) = 1.6a;o. (a) The incoming pulse spectrum 
\ex \ is located at a spectral region where |e| ~ 0.9 and therefore the interaction is outside of the ENZ regime, (b) For all the 
incoming pulses with amplitudes in the range 10“® < eo < 1.6 the absolute values of the normalized analytic signals |ef |/eo 
and |ef l/eo of the outgoing pulse are effectively identical, (c) For the pulse with amplitude eo = 1.2 the square absolute values 
|pf 1^ + IPz 1^ of the dimensionless analytic signal of the polarization within the slab are uniformly much smaller than 1. For the 
same pulse amplitude the dimensionless displacement field components spectra, \d^\ and |, coincide with the dimensionless 
electric field components spectra multiplied by permittivity absolute value, |eef^| and |eef |, within the slab. This results shows 
that pulse scattering outside of the ENZ regime is fully linear. 


nonlinear dependence of the scattering process on the in¬ 
put intensity saturates as expected from Eq. (HD. 

In Fig.2c the transition from the linear regime to the 
nonlinear one is pictorially shown in two different and 
equivalent ways: (i) the analysis of the polarization of 
the slab, and (ii) the spectral analysis of the field dy¬ 
namics within the slab. In the first row of Fig.2c we 
report the space-time profile of the polarization pro¬ 
duced within the slab by some of the pulses considered 
in Fig.2b, by plotting the corresponding square absolute 
values |pf(X,0,T)p -I- |pf(X,0,T)p of the dimension¬ 
less analytic signal of the polarization at the slab mid¬ 
dle plane Z = 0. The relevance of such quantity stems 
from the fact that it is comparable with |P|^ /which 
plays a key role in the linear-nonlinear transition of the 
wave-matter coupling (see Eq.([T|)) through the nonlinear 
resonant frequency 

. (4) 

For Co = 10“® the normalized square polarization is so 
small that We(P) — Wg. Accordingly, the linear regime 
(where the dielectric permittivity of Eq.(I5]) governs elec¬ 
tromagnetic propagation) holds. For Cq = 2 • 10“^ the 
normalized square polarization is small enough to allow 
a perturbative description of medium nonlinearity since 
a;g(P) ~ 1 — 3 |P|^/(2Pg^). In this case Eq.([T]) is well- 
known to yield a nonlinear optical regime characterized 
by a Kerr nonlinearity [b^. For cq = 1.2 the normalized 
square polarization is accordingly larger and does not 
generally allow a perturbative approximation of Wg(P), 
thus triggering the emergence of nonlinear saturation. 
In the second and third row of Fig.2c we plot the ab¬ 


solute values of the Fourier transforms d^{X,0,oj) and 
df {X, 0, oj) (green surfaces) of the components of the di¬ 
mensionless displacement field d = e -t- p at the slab mid¬ 
dle plane Z = 0. In the same rows we also plot the ab¬ 
solute values |e(a;)e,^(A, 0,a;)| and |e(a;)ef (A, 0, w)| (red 
surfaces) of the Fourier transforms of the electric fields 
components at the same plane Z = 0 multiplied by the 
linear permittivity of Eq. For eo = 10 ® it is evident 
that d'^ = ee^ so that the standard linear relation be¬ 
tween the displacement and the electric field holds. For 
eo = 2 • 10“^ and Cq = 1.2 the discrepancy between 
the displacement field and its linear counterpart is dra¬ 
matically evident thus restating that a marked nonlinear 
regime holds. 

The most remarkable nonlinear trait of the scattering 
at hand in the ENZ regime is that the normalized am¬ 
plitude of the outgoing pulse grows as the amplitude of 
the incident pulse is increased, as shown in Fig.2b. This 
phenomenology admits a simple explanation in terms of 
an effective nonlinear shift of the resonant frequency pro¬ 
duced by the emergence of the nonlinear regime. From 
Fig.2c it is evident that the higher the amplitude eo of 
the incoming pulse the stronger the polarization induced 
within the slab. From Eq.(|3]) this produces a local de¬ 
crease of the nonlinear resonant frequency Qe which can 
be equivalently interpreted as an effective drift of the 
pulse spectrum toward the spectral region where the real 
part of the dielectric permittivity is positive (see Fig.lb). 
Therefore as eo increases the slab becomes more trans¬ 
parent to the pulse and allows larger energy transmission. 
Since the pulse is quasi-monochromatic and its spectrum 
is centered at the zero-crossing-point, even a slight spec¬ 
tral drift produces a relatively large change of the mate¬ 
rial response, which explains the dramatic dependence of 
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FIG. 4: Mechanism supporting the full potential of the nonlinear wave-matter interaction in the ENZ regime. 

Spatial profiles of the dimensionless electric and polarization fields at the time T = 753 of the pulse with eo = 0.4 and 
tj = 1.6u!o (outside of the ENZ regime) and the pulse with eo = 4 • 10“^ and 4) = wo (in the ENZ regime). In the first case the 
electric field Cx is ’’large” only at the regions around the peaks of the wave and the time drift of such regions does not allow the 
polarization px to increase and to trigger the nonlinear wave-matter coupling. In the second case, in a physically large volume, 
the spatially slowly varying character of Cx yields the onset of the nonlinear regime. 


the slab transmissivity on eo reported in Fig.2b. 

It is worth noting that the interaction regime has a 
highly nonlinear character since a distinct pulse self- 
action occurs even if the slab is very thin, its thick¬ 
ness being comparable with the pulse carrier wavelength 
(L = 1.368Ao where Aq = 27rc/wo is the wavelength of the 
zero-crossing-point). In addition, we observe that the de¬ 
scribed nonlinear wave-matter interaction is not due to 
the enhancement of the longitudinal electric field compo¬ 
nent Cz Indeed, the transverse magnetic pulse 

incident from vacuum has a longitudinal component re¬ 
sulting from the finite size of the beam in the transverse 
direction (along the a;-axis) and, due to the field match¬ 
ing at the slab edge, the longitudinal component within 
the slab is roughly |e|“^ ~ 20 times larger than its vac¬ 
uum counterpart. However, from the second and third 
row of Fig.2c, it evident that |eea;| and jeczl are compara¬ 
ble within the slab so that \ex\ and |ez| are comparable as 
well. This demonstrates that in this case the longitudinal 
field component Cz does not play the usual leading role we 
have grown accustomed to, notwithstanding the fact that 
the longitudinal component is still enhanced. As a conse¬ 
quence, the same nonlinear mechanism is triggered with 
incident transverse electric pulses. Therefore, a distinct 
advantage one gains is that, by focusing more tightly an 
incident transverse magnetic pulse, both enhancement of 
the longitudinal field and the kind of nonlinear behav¬ 
ior we have just discussed can in principle combine for 
optimal results. 

Absence of nonlinear effects outside of the ENZ 
regime. The remarkable impact of matter nonlinearity 
on electromagnetic propagation in the ENZ regime is ren¬ 
dered even more peculiar and fascinating by the fact that 


no deviation from linear behavior occurs in the pulse scat¬ 
tering outside of the ENZ regime. This important fact 
is proved by a second set of simulations where we have 
considered various pulses identical to those discussed in 
Fig.2, except for their carrier frequencies, which we have 
set to w = I.Bwq. The results are reported in Fig.3, where 
we have plotted the same quantities considered in Fig.2. 
Figure 3a clearly shows that the incoming pulse spec¬ 
trum \e^\ is located in a spectral region where |e| ~ 0.9. 
Therefore, interaction does not occur in the ENZ regime, 
and the slab shows purely linear dielectric behavior. Fig¬ 
ure 3b shows that the pulse amplitude eo has no impact 
on the scattering since for all the considered pulses with 
10“® < eo < 1.6 the profiles of the absolute values of 
the normalized analytic signals |ef|/eo and |ef|/eo are 
all the same. Figure 3c reveals that the square absolute 
values -I- of the dimensionless analytic signal 
of the polarization within the slab, even for the relatively 
large pulse amplitude cq = 1.2, is uniformly much smaller 
than 1. In addition. Fig. 3c shows that, for the same 
large pulse amplitude eo = 1.2, the standard linear rela¬ 
tion = ee^ holds between the displacement and the 
electric field within the slab. As a result, pulse scat¬ 
tering does not reveal any nonlinear traits if the pulse 
carrier frequency is not close to the zero-crossing-point, 
no matter how large the pulse amplitude eo is made. We 
conclude that the slab governed by Eq. m does not play 
host to a nonlinear wave-matter interaction, unless the 
ENZ regime is exploited. 

Mechanism supporting the highly nonlinear 
wave-matter coupling in the ENZ regime. The 

discussed nonlinear wave-matter interaction that occurs 
only in the ENZ regime admits an explanation which is 
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closely related to the slowly varying character of the elec¬ 
tromagnetic field imposed by the very small permittivity. 
In Fig.4 we plot the spatial profiles of the dimensionless 
electric Cx and polarization px fields, at a specific time 
T = 753, for two of the pulses considered above, and hav¬ 
ing the same amplitude eg = 4-10“^ and different carrier 
frequencies w = 1.6wo and Co = loq. The first pulse (left 
column of Fig.4) is tuned outside of the ENZ regime. 
Its scattering is purely linear, whereas the second pulse 
(right column of Fig.4) with (D = wq is in the ENZ regime 
and it has been shown to display a marked nonlinear dy¬ 
namics. Consider the first pulse with Co = I.Gwq and note 
that, at the beginning of the scattering interaction, pulse 
propagation is evidently linear since the initial electric 
field is small. Therefore the profile of the pulse electric 
field has a considerable number of nodes, the regions of 
low electric field around such nodes drift in time due to 
propagation, the amplitude |A(r,t)| of the electric field 
E(r, t) = Re [A(r, driving the polarization at 

each point rapidly varies between zero and it maximum 
value. Since in the linear regime Eq. © yields 

V{t)=egf dt'x{t-t'W), (5) 

where x(t) is the Lorentz susceptibility, we conclude that 
the rapid variation of the electric field amplitude forbids 
the polarization from locally increasing and from driving 
the scattering process out of the linear regime. In the 
case of the second pulse with Co = uog the initial linear 
dynamics forces the electric field to be spatially slowly 
varying within the bulk since it is in the ENZ regime. As 
a consequence, very few nodes appear in the electric field 
profile and correspondingly the amplitude |A(r, t) \ of the 
electric field E(r, t) = Re [A(r, has a time vari¬ 

ation scale much slower than the first pulse, occuring over 
a physically large portion of the bulk. Hence the polar¬ 
ization is efficiently pumped by the electric field (much 
more than the pulse tuned outside of the ENZ regime) so 
that it correspondingly increases to the point of trigger¬ 
ing the onset of the nonlinear wave-matter coupling we 
have discussed above. 


II. DISCUSSION 

In conclusion we have shown that a novel and highly 
nonlinear wave-matter coupling occurs if the medium’s 
linear permittivity is very small over the entire electro¬ 
magnetic field bandwidth. The realization of the full po¬ 
tential offered by this nonlinear matter response is due to 
the spatially slowly varying character of the electromag¬ 
netic field in the ENZ regime. In fact, the ENZ regime 
produces more efficient coupling between the electric field 
and the medium polarization field by enlarging the effec¬ 
tive portion of the bulk hosting the nonlinear interaction. 
Unlike most strategies proposed to date in literature in 
order to achieve a highly nonlinear response, the nonlin¬ 


ear regime we have discussed here is not triggered by ei¬ 
ther cavity effects or by large nonlinear coefficients. The 
strategy we have outlined is simpler and more suitable to 
be exploited in a number a different configurations since 
it requires only that the field’s main frequency coincides 
with a zero-crossing-point of the real part of the dielec¬ 
tric permittivity, preferably at the crossing point where 
absorption may be neglected. In addition, in view of such 
flexibility, the nonlinear coupling we have discussed may 
be triggered even in the presence of other mechanisms 
that are known to enhance matter-wave coupling, thus 
likely to produce additional hitherto unknown, highly 
nonlinear effects. For example, the enhancement of the 
field component normal to the surface between vacuum 
and the ENZ material may be triggered in the scattering 
experiment considered in this paper simply by confining 
more tightly the incoming pulse along the transverse di¬ 
rection. The simplicity, generality and flexibility proper¬ 
ties of the discussed highly nonlinear regime make it an 
ideal platform for conceiving a number of applications 
and a novel generation of ultra-compact and fast devices 
for manipulating light. 


III. METHODS 


Full -wave simulations. Electromagnetic propaga¬ 
tion within the slab is described by Maxwell equations 
coupled to Eq. o for the polarization so that, by us¬ 
ing the dimensionless coordinates R = Wef/c, T = lOet 
and dimensionless fields e = egE/Ps, h = H/(cPs), 
p = P/Ps, the pulse scattering is described by the set 
of equations 


Vr X e 


dp 


Vr X h 



(9h 

de dp 
(es-l)e (6) 


within the slab, by the first two equations with p = 0 
outside of the slab and matching conditions of the tan¬ 
gential electric and magnetic field components at the slab 
surfaces. 


The pulse scattering by the slab both in and outside 
of the ENZ regime was simulated by solving the above 
equations through a home-made finite-difference time- 
domain code suitable for dealing with transverse mag¬ 
netic pulses e{X,Z,T) = ex{X, Z,T)i + ez{X, Z,T)'k, 
\i{X,Z,T) = hy{X, Z,T)}. At the edges of the com¬ 
putational domain perpendicular to the Z axis scatter¬ 
ing boundary condition where adopted whereas at those 
perpendicular to the X axis the vanishing of all the field 
components was imposed. 
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